Field-induced magnetism in the high-T c superconductor Ca 2 _ :r Na. r Cu02Cl2 

with flat Cu0 2 planes 



K. Ohiahi/'B I. Yamada, 2 ^ A. Koda, 1 - 3 S.R. Saha, 1 '!! R- Kadono, 1 * 3 
W. Higemoto, 4 K.M. Kojima, 5 M. Azuma, 2 and M. Takano 2 -§ 

institute of Materials Structure Science, High Energy Accelerator 
Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 
2 Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan 
'Department of Materials Structure Science, The Graduate University 
for Advanced Studies (SOKENDAI), Tsukuba, Ibaraki 305-0801, Japan 
'Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan 
5 Department of Physics, Graduate School of Science, 
University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
(Dated: October 30, 2009) 

The internal magnetic field distribution in a mixed state of a cuprate superconductor, 
C^-aiNa^CuChCk (T c — 28.5 K, near the optimal doping), was measured by muon spin rota- 
tion (/iSR) technique up to 60 kOe. The piSR linewidth A(B) which exhibits excess broadening 
at higher fields (B > 5 kOe) due to field-induced magnetism (FIM), is described by a relation, 
A(B) oc \TS. This suggests that the orbital current and associated quasiparticle excitation plays 
predominant roles in stabilizing the quasistatic correlation. Moreover, a slowing down of the vortex 
fluctuation sets in well above T c as inferred from the trace of FIM observed up to ~ 80 K, and 
develops continuously without a singularity at T c as the temperature decreases. 

PACS numbers: 74.25.Ha, 74.72.-h, 76.75.-K 



The 1/8 anomaly, or the suppression of superconduc- 
tivity at a specific hole carrier concentration, p ~ 1/8, 
has been understood in terms of the so-called "stripe" 
model, where a quasistatic charge and spin stripe order 
develops in place of superconductivity^. While a sim- 
ilar situation is realized over a wide range of p ~ x 
in La 2 _ £! ;_yNdj / Sr a ;Cu04 (RefQ]), a spatially modulated 
dynamical spin correlation is found in L^-xSr^CuC^ 
(LSCO) that becomes static only near p ~ 1/8 (RefH. 
It is speculated that the dynamical stripes tend to be 
static at p ~ 1/8 by a synchronization of the periodic- 
ity between lattice and charge/spin modulations, as the 
buckling of CuC>2 planes in the tetragonal structure ob- 
served at low temperatures (LTT phase) seems to serve as 
a source of "pinning" , leading to a stronger 1 /8 anomaly. 

Recent neutron scattering experiments in under- 
doped and optimally doped LSCO^^ and underdoped 
La2-a:Ba K Cu04 (LBCO)£ have revealed that the qua- 
sistatic stripe order is enhanced at around p ~ 1/8 by 
applying a moderate magnetic field of a few tesla paral- 
lel to the c axis, as inferred from the enhanced intensity of 
incommensurate magnetic Bragg peaks. Moreover, such 
a field-induced static magnetic order has turned out to 
exhibit a three-dimensional correlation, which survives 
into the superconducting state, in LSCO (p = 0.10), sug- 
gesting a predominant role of the orbital current in field- 
induced magnetism (FIM) 7 -. 

The occurrence of FIM in high-T c cuprates, which in- 
cludes static and/or dynamical magnetic correlations, 
has also been suggested from muon spin rotation 
(/iSR) measurements in underdoped and optimally doped 
LSCOSAmi, LBCC£, La 2 _ x _ y Euj / Sr x Cu04 (LESCO)2, 
and YBa2Cu307_5 (YBCO) 12 . There is a common fea- 



ture that the muon spin relaxation is observed well above 
T c in these cuprates, suggesting that the FIM develops 
well above T c . Recently, the FIM has been observed even 
in overdoped LSC 12 i 13 ' 14 . The field-induced static mag- 
netic order observed by neutron scattering experiments, 
on the other hand, occurs only below T ~ T c in under- 
doped and optimally doped samples of LSCO. Here, the 
field-induced static magnetic order is presumed to origi- 
nate from quasistatic moments around the vortices, and 
such a magnetic order is not observed for overdoped sam- 
ples. A lack of experimental evidence from neutrons for 
static magnetic order persisting above T c led to the ar- 
gument that the quasistatic FIM observed above T c by 
/iSR could be due to uncorrected static spins. 

Regardless of such a correlation with T c , the field- 
induced static magnetic order has a common feature that 
it is observed in those high-T c cuprates that exhibit a 
buckling of the Cu02 planes. From this view point, 
Ca 2 _ 2; Na a ;Cu02Cl2 (Na-CCOC) is a good candidate for 
studying the relation between the FIM and the buckling 
of Cu02 planes, because Na-CCOC has flat Cu02 planes 
without buckling, and therefore it would serve as a stage 
for testing the FIM as an intrinsic character of Cu02 
planes. 

In this paper, we demonstrate the occurrence of FIM in 
Na-CCOC probed by //SR. While the crystal structure of 
Na-CCOC is isostructural with that of LSCO, it consists 
of flat Cu02 planes, owing to the substitution of apical 
oxygen with chlorine. While the /iSR linewidth A in the 
mixed state decreases with increasing field at low mag- 
netic induction, as expected for a normal flux line lattice 
(FLL), it exhibits a turnover at around B ~ 5 kOe, and 
increases in proportion to \fB, strongly suggesting that 
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FIG. 1: (Color Online) Fast Fourier transform of TF-^SR 
spectra observed at 40 K (solid curve, H = 0.5 kOe) and 2 K 
(dashed curve, 0.5 kOe; dotted curve, 60 kOe). 



the depolarization is due to quasistatic magnetism asso- 
ciated with magnetic vortices whose quasiparticle density 
is proportional to y/B in d-wave superconductor o 15 ' 16 . It 
is inferred from the muon Knight shift measurement that 
local spins of 0.15^b/Cu are responsible for the observed 
effect. Moreover, an enhancement of A is observed up 
to 80 K under a field of 60 kOe, which is well above T c 
(=28.5 K). This reveals the possibility that the fluctu- 
ation of random vortices persists over a wide range of 
temperature far exceeding T c , which is in intriguing ac- 
cord with the Nernst effec t 17 > 18 . 

The single-crystalline specimen of Na-CCOC used in 
this study had been grown by a flux method under high 
pressure to yield slab samples with the c axis normal 
to their plan o 19 i 20 . The superconducting transition tem- 
perature determined from susceptibility measurement is 
28.5 K (see the inset of Fig. 2]), which corresponds to 
that of optimally doped compounds according to the 
phase diagra m 21 i 22 . The samples were encapsulated in 
a polymide tape in a glove box to prevent the depletion 
of sodium and a subsequent deterioration. Transverse 
field (TF) /iSR measurements were performed on the M15 
beamline of TRIUMF, Vancouver, Canada. Muons with 
a momentum of 29 MeV/c were injected into the samples 
with their polarization rotated perpendicular to the beam 
momentum, so that the external field H might be applied 
along the incoming beam axis (to minimize the distur- 
bance to muon trajectory), which was parallel to the crys- 
talline c axis. An experimental setup with a high time 
resolution was employed to measure TF-/iSR time spec- 
tra up to 60 kOe. The corresponding dc-susceptibility at 
60 kOe was measured by a SQUID magnetometer. 

Figure Q] shows the fast Fourier transform (FFT) of 
TF-/iSR time spectra, which represent the density distri- 
bution of internal magnetic field B at the muon site (with 
additional broadening due to random local fields from 
nuclear magnetic moments and that coming from a finite 
time window of FFT). In the normal state above T c (solid 
curve), the spectrum has a peak at B = H = 0.5 kOe, 



which undergoes a shift to a lower frequency upon cool- 
ing to 2 K with an increased linewidth (dashed curve). 
This behavior is perfectly in line with the FLL formation 
in the superconducting state. Meanwhile, the spectrum 
exhibits further broadening with increasing field to 60 
kOe (dotted curve), which is opposite to the predicted 
tendency of decreasing linewidth with increasing field in 
the conventional type II superconductors^. The data 
are analyzed by curve fits in the time domain using a 
phcnomcnological stretched exponential function, 



AP(t) = Aexp -{Atf exp(zw At t 



(1) 



where A is the positron decay asymmetry, A is the depo- 
larization rate (= linewidth in the frequency domain), (3 
is the power of the exponent, u>^ = j^B with 7^ being 
the muon gyromagnetic ratio (= 27rx 13.553 MHz/kOe), 
and <fi is the initial phase of precession. The muon Knight 
shift, K, is then defined as, 
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where ojq — Jfj,H, is the shift due to the local 
spin susceptibility Xspin, K v is that due to the orbital 
current in the FLL state, and K dcm is the correction 
term consisting of demagnetization and Lorentz field 
[= 47r(l/3 — AQpmolXmolj where N ~ 1 is the demag- 
netization factor, p mo i = 0.01480 mol/cm 3 is the molar 
density, and x mo i is the molar susceptibility]. 

Figure [2] shows the field dependence of the relaxation 
rate at 40 K and 2 K, and that of the power of the ex- 
ponent at 2 K. The relaxation rate, a (= A^=2), at 40 K 
exhibits the iiz-linear behavior, consistent with previous 
results in LSCC&i^ L BCC£, LESCQ2, and YBCOi 2 -. 
Note that this linear behavior of a can be understood by 
assuming the existence of slowly fluctuating (staggered) 
random magnetic moments induced inside the vortices. 
As mentioned later, it is described that these fluctuating 
moments appear with the fluctuation of vortices, because 
the number of vortices corresponding to the fluctuating 
moments is proportional to the applied field. Therefore, 
this iJ-linear behavior of a is considered to be slowly 
fluctuating random magnetic moments. On the other 
hand, while A at 2 K decreases with increasing field at 
lower fields, it exhibits a turnover around H = 5 kOe 
and an increase represented by the relation A cx VH 
(where H ~ B). The lineshape shows a change from 
that of the Gaussian ((3 = 2) to single exponential de- 
cay ((3 = 1) with increasing field, as shown Fig. [2jb). 
The contribution of FLL at low fields, cr v , is extracted by 
subtracting a n from the total linewidth, <r, in quadrature, 
(T 2 = cr 2 — fi 2 , with a n being the depolarization rate in the 
normal state (see below). The behavior observed at lower 
fields is consistent with the predicted //-dependence 2 ^, 
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(1 -h)y/l + 3.9(1 -h) 2 , (3) 



where $0 is the magnetic flux quantum, A t(T, h) is the 
effective inplane London penetration depth, which can 
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FIG. 2: (Color Online) Magnetic field dependence of (a) the 
relaxation rate A and a at 2 K and 40 K, respectively, and 
(b) the power of exponent /3 at 2 K. Inset shows the \f~H 
dependence of A above H = 10 kOe. Solid curve represents a 
fit by the relation A oc vH. 



be expressed as A a &(T, 0) ■ (1 — rjh) with -q > for line 
nodes^, and h is the field normalized by the upper crit- 
ical field (h = H / H c2 ). Meanwhile, the value of A above 
~5 kOe far exceeds that observed at lower fields. It is 
unlikely that such an enhancement at a higher field is 
induced by flux pinning or other extrinsic artifacts, and 
thus can be uniquely attributed to FIM that is spatially 
inhomogeneous (as inferred from the single exponential- 
like lineshape). More interestingly, A is excellently re- 
produced by the relation A(B) oc \f~B over the relevant 
field range. The inset of Fig. [2](b) shows A, which is 
plotted against \f~B to see the linearity. Note that this 
behavior is not observed in other high-T c cuprates that 
do not have flat Cu0 2 planes&±&±i>±^>±l It is estab- 
lished that quasiparticle excitation in the mixed state 
of d-wave superconductors is extended along the (it, n) 
directions, leading to a non- linear field dependence of 
the quasiparticle density that is well approximated by 
such -\/B-dependenc o 15 ' 16 . Here, we would like to con- 
sider the relation between the muon spin relaxation and 
the quasiparticle excitations in d-w&ve superconductors. 
In the case of the usual type II superconductors with 
(i-wave symmetry, which do not show any FIM, muon 
spin relaxation rate decreases with increasing H . This 
is because the induced quasiparticle excitations around 
the nodal region due to pair breaking makes A increase. 
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FIG. 3: (Color Online) Temperature dependence of the re- 
laxation rate (a) a and ov (inset) at H = 0.5 kOe, (b) A 
at H = 60 kOe, and (c) the total muon Knight shift K at 
H — 60 kOe. The solid curve in the inset of (a) is a fit with 
CTv(T) = <t(0)[1 - (T/T c )"]. Those in (b) and (c) are the 
guides for the eyes. 



According to a relation of Eq. ([3]) , a v decreases with in- 
creasing A. On the other hand, in the case of super- 
conductors, which show FIM, muon spin relaxation rate 
increases because relaxation due to magnetic moments 
is dominant compared with that due to FLL (increase 
of A)£&±&±ii±^i±l It is suggested that the induced 
quasiparticle excitations have quasistatic magnetic mo- 
ments to describe the increase of A(B). Consequently, 
we can assume that the observed field dependence of A 
is one type of evidence for the d-wave superconductor, 
because quasiparticle excitations in d-wave superconduc- 
tors shows a -\/B-dependence. 

As shown in Fig. [3^a), while the linewidth, a (— A ( 3 = 2 
observed at 0.5 kOe), is mostly independent of tempera- 
ture above T c with a mean value of a n = 0.144(3) /J.s _1 , it 
increases with decreasing temperature below T c as FLL is 
formed. The value of a n is in perfect agreement with that 
of random local fields from nuclear moments estimated 
from a previous /xSR experiment 2 —. Fits by a power law, 
cr v (T) = cr v (0) [1 - (T/T c )"], with T c as a free parame- 
ter yields cr v (0) = 0.519(5) ^ , T c = 27.7(2) K, and 
n = 3.2(1) [Fig. G2a), inset]. Using Eq. © for h < 1, 
the magnetic penetration depth, A a fc(0, 0), extrapolated 
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FIG. 4: (Color Online) (a) Temperature dependence of the 
muon Knight shift, (the spin part, red circle), and sus- 
ceptibility (blue triangle) at H = 60 kOe. Inset shows the 
temperature dependence of susceptibility observed at H = 10 
Oe. For the origin of negative offset, see text, (b) versus x 
plot. The turnover below ~10 K is attributed to unidentified 
impurities. 



to T = K is evaluated to be 382(4) nm. These re- 
sults are quantitatively consistent with earlier literature 
on Na-CCOC (with the sodium concentration x = 0.18, 
and H = 2 kOe), except that no anomalous increase of 
a was observed below ^5 K 25 . 

The temperature dependence of A observed at H = 
60 kOe is shown in Fig.^b). It increases below 80 K (3> 
T c ), and continues to increase down to the lowest tem- 
perature without any saturation. This behavior is simi- 
lar to those observed in LSCO at a high field2i±2 i ilil2il£. 
On the other hand, the muon Knight shift, K, decreases 
with decreasing temperature below T c [Fig. G3c)], indi- 
cating that both the FIM and superconductivity coexist 
below T c . The behavior of K is in line with that of the 
magnetic susceptibility [%, measured after cooling under 
a field of 60 kOe, shown in Fig. Ufa)], except that below 
~10 K (see below), although the contributions of spin 
and orbital parts are not separated in these quantities. 
The spin part of the muon Knight shift, K^, is extracted 
by Eq. © for T > T c (where K v = 0) and mapped 
into the K^-x plot [Fig. [Hb)]. All of the data points 
fall on a straight line without any sign of kink for those 
obtained below ^80 K (smaller values of K^-x) , indicat- 
ing that FIM is indeed carried by the paramagnetic Cu 
spins. The muon hyperfine parameter, A^, is deduced as 



TABLE I: Summary of parameters for muon sites identified in 
Na-CCOC, where is the fractional yield, rc u is the distance 
between muon, A% z is the dipole tensor calculated for the 
muon coordinates shown in the preceding column. The crys- 
tal structure of Na-CCOC is K2NiF4-type structure (space 
group: IA/mmm) with lattice constants a — 3.8687A and 
c = 15.0485A. 
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a gradient, dK^/dx, in the normal state according to, 
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which yields A^ — 696(2) Oe/MB, where K (~ 0) is a 
T-independent Fermi contact coupling with the conduc- 
tion electrons, 2Va is the Avogadro's number and Mb is 
the Bohr magneton. The mean value |pc u | f° r the mo- 
ment size of the local copper spins is then obtained by 
comparing A^ with that determined by the z-component 
of the dipole tensor, 
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where is the relative population of the fc-th muon site 
(there are three of them in Na-CCOC 2 ^, see below), 
is the distance between muon at the k-th site and the 
i-th Cu ions (with Zik being their z-component). Com- 
bining A^ with other information on the local magnetic 
fields at the muon sites previously obtained for antiferro- 
magnetically ordered phase of the parent compounds^, 
we reexamined the muon sites to reproduce all of the 
quantities concerning the local fields. The result is sum- 
marized in Table [H where the /ii site turns out to be 
slightly off the previously assigned position as inferred 
from the sign of A^. Accordingly, we have A dip = 
^2k a kA% z = 4788 Oe/MB, and the effective moment size 
|7*c3ul = 0.1454(4)mb from Eq. ©. 

It is noteworthy in Fig. [4^a) that x exhibits an upturn 
below ~10 K, which might hint the occurrence of mag- 
netic order. However, no such behavior is observed in the 
muon Knight shift K shown in Fig.[3^c); it must be noted 
that (the spin part) is deduced using Eq. ([2|) where 
X comes in via the demagnetization term (Adcm oc x)i so 
that the temperature dependence of is mostly parallel 
with that of x- We attribute this upturn in x to unidenti- 
fied impurities at this stage, considering the absence of a 
corresponding anomaly in A^ anticipated for quasistatic 
magnetic order. The ambiguity on the absolute values of 
X makes it difficult to discuss whether or not the observed 
change of sign in below ~20 K is meaningful. How- 
ever, such a behavior is readily explained by considering 



5 



the contribution of K v (< 0) which is not discernible 
from other contributions in the present experiment. 

Since Na-CCOC does not exhibit a quasistatic stripe 
correlation under a zero external field over the relevant 
hole concentration 2 ^, the occurrence of FIM can not 
be explained by a residual effect of the 1/8 anomaly, 
as it might have been in LSC02&. A recent reve- 
lation of a weak 1/8 anomaly observed in Na-CCOC 
upon Zn substitution for Cu^ 7 - strongly suggests that 
the situation is similar to that observed in YBCO 28 
and Bi 2 Sr 2 Cai_ 2; Y ;l; Cu20 8+ ^. This, together with the 
present result, supports the presumption that the insta- 
bility of dynamical stripe correlation against local sup- 
pression of superconducting order parameter (i.e., by the 
Zn impurity or flux lines) is an intrinsic property of Cu02 
planes, common at least in Na-CCOC, LSCO and LBCO 
regardless of the buckling of the Cu02 planes. 

According to the reported result on the Nernst 
effec t 17 ' 18 , it is suggested that strong fluctuation of super- 
conductivity and associated flux lines persists well above 
T c . This naturally leads to a possibility that the FIM 
observed above T c may originate from the quasistatic 
stripe correlation segregated around the fluctuating vor- 
tices. Considering the result of field dependence of A at 
40 K, which behaves 77-linear, this scenario is apparently 
against the expected field dependence in the FLL stated. 
However, it may be argued that the V^B-dependence is 
expected only when the quasistatic FLL is established 15 . 
Moreover, this scenario is also supported by thermal con- 
ductivity measurements in LSCO, where the suppression 
of thermal conductivity observed above T c is attributed 
to the development of the quasistatic stripe order asso- 
ciated with fluctuating vortices^. Thus, the appearance 
of FIM above T c can be explained in terms of the insta- 



bility of dynamical stripes upon the local suppression of 
the order parameter. 

Finally, it would be worth mentioning that the sup- 
pression of the FIM near p ~ 1/8 observed in YBCO 12 
has little or no direct relevance to our argument. While 
the relaxation rate is suppressed near p ~ 1/8 compared 
to other hole concentrations, they indeed observe that the 
FIM sets in at around 210 K at p ~ 1/8, which is much 
higher than T^. Thus, YBCO seems to share the fea- 
ture of FIM that develops well above T c (while the onset 
temperature for the FIM might exhibit the p-dependence 
similar to T c characterized by a dip near p ~ 1/8). 

In conclusion, our ^SR measurements on optimally 
doped Na-CCOC demonstrate the occurrence of FIM 
above ~5 kOe, which is highly inhomogeneous and co- 
exists microscopically with superconductivity. Muon 
Knight shift measurements indicate that the magnetism 
is carried by local copper spins with an average moment 
size of 0.15^b- The field-dependence of the muon spin 
depolarization rate suggests that the magnetism comes 
from the quasistatic stripe correlation around the vortex 
cores (with ci-wave paring). This is also in line with the 
interpretation that the observed FIM up to 80 K (well 
above T c ) comes from the strong fluctuation of vortices 
suggested by the Nernst effect. 
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